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We show that relativistic mean fields theories with scalar, S, and vector, V, quadratic 
radial potentials can generate a harmonic oscillator with exact pseudospin symmetry 
and positive energy bound states when S = —V. The eigenenergies are quite different 
from those of the non-relativistic harmonic oscillator. We also discuss a mechanism for 
perturbatively breaking this symmetry by introducing a tensor potential. Our results 
shed light into the intrinsic relativistic nature of the pseudospin symmetry, which might 
be important in high density systems such as neutron stars. 
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Recently, there has been much interest in understanding the nuclear pseudospin 
symmetry in terms of relativistic dynamics . However, in the presence of nuclear 
scalar and vector mean-field potentials this symmetry is broken non-perturbatively 
as it has been discussed in recent articles 2,a . In this work we will show that for 
quadratic radial scalar, S, and vector, V, potentials with S = — V is possible to 
break the pseudospin symmetry by introducing a tensor potential that preserves 
the form of the harmonic oscillator central potential but generates a pseudospin 
orbit term. 

1. Harmonic oscillator and exact pseudospin symmetry 

In relativistic mean field theories with scalar and vector potentials the Dirac 
Hamiltonian for a fermion (nucleoli) of mass m is 

ff D = a.p + /3m + i(l + /3)S + i(l-/3)A , (1) 
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where a and (3 are the usual Dirac matrices. We have introduced the "sum" and 
the "difference" potentials defined by £ = V + S and A = V - S. When E = or 
A = 0, due to the matrix structure of ^(1 ± j3), the Dirac Hamiltonian is invariant 
under a SU(2) symmetry . This is a general feature which does not depend on the 
particular form of S and V. In this paper we will consider only the case £ = 
(S = —V) and a harmonic oscillator potential with angular frequency u\ for the A 
potential, i.e., A = \rrvJ 2 r 2 . 

A more generalized form of the relativistic harmonic oscillator where the case 
A = is also discussed has been presented recently 4 . Following the details of that 
paper it is easy to prove that, for the case £ = 0, the lower component of the Dirac 
spinor satisfies the differential equation 



d _ 1(1 + 1) _ rn(£-rn)_ 
dr^ r z 2 



/ [ I — !— II 111 I I i I I 

f K (r)=0. (2) 



As shown in , I is not the non-relativistic angular momentum I but rather 
the pseudo-orbital angular momentum of the lower component of the Dirac spinor. 
They are related by £ = £ — k/\k\, where the quantum number k determines whether 
spins are parallel or antiparallel. Pseudospin symmetry is exact when doublets with 
j = I + s are degenerate, where s = s is the pseudospin quantum number. 

The eigenenergies are given by 



(S + m)\ ?—^ =wi ( 2n + l + \ ) , (n = 0, 1, 2, . . .) . (3) 



2m V 2 

where n is the radial quantum number related to the nodes of the lower component 
f K (r) . From Eq. (3) one concludes that the real solutions must have positive binding 
energy E = £ — m. The second order equation (2), which only depends on /, and 
also the eigenenergy equation (3), show that there is no pseudospin-orbit coupling. 
Therefore the states with same (n, I), but with j = 1 + 1/2 and j = I — 1/2 are 
degenerate: they can be seen as pseudospin doublets as discussed in 1 . Thus, when 
A is a harmonic oscillator potential and S = 0, pseudospin symmetry is exact and 
there are only positive-energy bound states, i.e., A acts as a binding potential. This 
is an interesting result in view of the fact that the pseudospin symmetry obtained 
in the limit S(r) ^0 cannot be realized for nuclear vector and scalar mean-fields 
which go to zero as r — > oo. In fact, in that case, X acts as a binding negative 
central potential well and therefore no bound states may exist when S = 1,2 . 
The spectrum of single particle states for the case E = 0, oj\ = 2 and m = 10 
is shown in Fig. 1 (a) using the quantum numbers of the upper component, that 
can be related to the non-relativistic quantum numbers. In Fig. 1 (b) we classify 
the same energy levels by the quantum numbers of the lower components f K . The 
comparison between these two figures exhibits the pseudospin symmetry and its 
quantum numbers. For example, the doublets [lsi/2 — 0d 3 / 2 ] and [lps/2 ~ O/5/2], 
which have the same pseudo angular momentum I and the same n, are, in the new 
notation, [Opi/2 — OP3/2] and [0^3/2 — CW5/2], respectively. Therefore, the harmonic 
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Figure 1: a) Single particle energies for the case S = with the quantum numbers 
(nlj) and b) (nZ j = j). The parameters are uo\ = 2 and m = 10. 



oscillator with £ = and A = ^muj\r provides an example of exact pseudospin 
symmetry. Fig. 1 (a) also shows the (2n + V) degeneracy, which means that not only 
the states with same n, I are degenerate (pseudospin partners) but also, for example, 
(n — 1, 1 + 2) or (n + 1, 1 — 2) have the same energy. 

The non- relativistic limit of the eigenvalue equation in Eq.(3) is reached when 
lo\ m, which, in turn, means that E = £ — m, <C 1, giving 



E 



1 

2^ 



2n + I + 



col 



(4) 



This equation, valid for T, = 0, shows that the energy is of second order in the 
ratio LJi/m, meaning that the energy is zero up to first order in LO\jm. We can 
interpret this fact by saying that, up to this order, there is no non- relativistic limit 
for S = and therefore the theory is intrinsically relativistic and so is the pseudospin 
symmetry. 

2. Perturbative breaking of the pseudospin by a tensor potential 

When X 7^ (S ^ — V) a pseudospin orbit term shows up in the equation for 
the lower component given by 

S * (5) 



£-£(r) r' 

where the prime denote the derivative with respect to r and k is the quantum 
number referred to before. So when S' ^ the denominator, E — X(r), can become 
very small near the pole and the contribution of this term can be quite significant, 
as shown in 5 . This manifests the non-perturbative breaking of this symmetry when 
S 7^ 0, as it is the case for the nucleus 2 . However, we may overcome this problem 
if we keep £' = £ = but a tensor potential i(3a • fU is introduced in the Dirac 
equation. In this case a new term appears in the spin-orbit term, —2nU/r, as shown 
in 4 . Because of the product a • r, the pole structure in the denominator discussed 
before is absent. 
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Figure 2: Energy levels for £ = with u)\ = 2 and m = 10 as a function of u 2 - 
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Figure 3: Energy spectrum for S = with u>i = 2, UJ2 = 1 and m = 10. 



If we choose a linear form for U (r) = mc^2 r, we have an additional quadratic 
central potential with a frequency w 2 in Eq. (2), so that we still have a global 
harmonic oscillator central potential. Moreover, we generate a pseudospin orbit 
term -2crau2. Thus, with this new term, the I degeneracy (pseudospin symmetry) 
can be broken perturbatively in the sense that can be made so small as the u>2 
frequency. We show in Fig. 2 how this symmetry is broken for the p states, for 
positive L02, and present the eigenencrgies in the spectroscopic notation nlj , where 
states with I = 0, 1, 2, . . . are denoted by s ,p , d , . . . respectively. Interestingly, from 
Fig. 3 we see that with the breaking tensor potential there is a value of u>2 for which 
Opi/2 (l s i/2): a state with I = rather than I — 0, becomes the ground state. 

Finally, the relativistic nature of the pseudospin symmetry with harmonic os- 
cillator potentials shown here suggests that this symmetry may be more important 
in ultra-relativistic systems. Thus, this symmetry may play a role in high density 
matter, such as the matter inside compact stars. 
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